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The strength of frustration and quantum fluctuations in LiVCuO i 
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Abstract -We present an empirical and microscopical analysis of the main in-chain exchange 
constants of the edge-shared frustratrated chain cuprate LiVCu04= LiCuVCU with a ferromag- 
netic nearest neighbour coupling which clearly exceeds the antiferromagnetic (AFM) next-nearest 
neighbour exchange Ji . The measured saturation field is significantly affected by a weak 3D AFM 
interchain coupling leaving room for a possible Bose-Einstein condensation for several T below. 
The obtained exchange parameters are in agreement with the results for a realistic five-band ex- 
tended Hubbard Cu 3d O 2p model, LSDA+17 predictions as well as with inelastic neutron and 
magnetization data. The single chain frustration rate a = J2/IJ1I ~ 0.75, including all error 
bars, is definitely smaller than 1.0 which correspond to strongly coupled interpenetrating AFM 
Heisenberg chains in contrast with opposite statements in the literature. A proper account of 
strong quantum fluctuations and frustration is necessary for a correct assignment of the exchange 
integrals, which cannot be achieved by a simple renormalization of J2 from spin-wave theory. 



1. INTRODUCTION. - LiVCu0 4 = LiCuV0 4 is 
one of the first [1,2] and rather frequently studied spin- 
chain compounds among edge-shared cuprates [3-7] . Re- 
cently it became especially interesting due to the observa- 
tion of multiferroicity [8-10] and due to a possible realiza- 
tion of quantum spin nematics and related Bose-Einstein 
condensation of two-magnon bound states in high mag- 
netic fields [11-13]. Both phenomena are still poorly un- 
derstood and a precise knowledge of the main exchange 
interactions is of key importance to attack such complex 
problems in a realistic way. Unfortunately, there is no con- 
sensus about the magnitude of these couplings and in par- 
ticular on the value of the in-chain frustration parameter 
a = J2/ 1 J\ I , with J \ < as the ferromagnetic (FM) near- 
est neighbor (NN) and J2 the antiferromagnetic (AFM) 
next-nearest neighbor (NNN)-coupling in chain direction 
b (see Fig. 1). So far, in various studies 0.5 < a < 2.2 and 
even above 5.5 have been predicted/reported. [3,5,14-16] 



Keeping in mind the weak interchain coupling, the single- 
chain can be viewed also as two interacting and interpen- 



t. (b) 




( a ) Corresponding author, E-mail: s.l. drechslerOifw-dresden.de 



Figure 1: (Color) (a): the crystallographic structure of 
LiVCuCU comprises two AFM coupled CuC>2 spin-chains per 
unit cell running along the 6-axis (orange • - Cu 2+ , red • - 
O 2- , bright blue - Li + ). (b): the main in- and inter-chain 
exchange paths, Ji, J2, and Js,J4, J5 marked by solid red arcs, 
broken red line, blue and green arcs , respectively. 
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Figure 2: (Color) Left: Magnetic dynamical structure factor S(k,uj) from exact diagonalizations for periodic chains and L — 28 
sites for the parameter sets proposed in Ref. [5] (a), our recent findings [16] (b) and that by Koo et al. [15] for [7 = 5 eV (see Tab. 
1) (c). Right: experimental asymmetric dispersion along the chain direction in between fl+ and Q,- (see Fig. 2(h) and Eq. (2)) 
and INS-intensities from Ref. [5] (d,e), calculated S(k,u>) normalized to the static structure factor S(k) (i.e. the w-integrated 
dynamical one) for a single chain with our parameters (f) and that of Ref. [5] (g), and the ratio of the peak positions at the 
first and second maximum of S(k,to) for the transferred momenta near fc=l/4 an 3/4, respectively, (h) approximately given by 
the red curves in (a-c). In all DMRG- calculations for chains with open boundary conditions, L = 96 sites, and a Lorentzian 
broadening at half-width F = 0.3 meV have been employed. 



etrating simple AFM Heisenberg chains (AHC) or equiv- 
alently as a single zigzag ladder. Then, one is left with a 
weak (a > 1) or a strong coupling scenario in the op- 
posite case < a < 1 because 1/a provides a direct 
measure of the FM "interchain" -coupling within a zigzag- 
ladder or between interpenetrating AHC. Here, we give 
a more comprehensive presentation of arguments against 



weak coupling scenarios than in our short Comment [16] 
on Ref. [5] and rebut in more detail also arguments put 
forward in Ref. [6]. In addition, the striking discrepancy 
of a very recent parameter set obtained in Ref. [15] with 
available experimental data, including also Rcfs. [5,6] will 
be shown, too. Finally, the reasons of the incorrect pa- 
rameter assignment given in Refs. [3, 5, 15] are explained 
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Table 1: Exchange integrals (in meV) as extracted from INS 
data using bare spin wave theory (BSWT), ad hoc renormalized 
spin wave theory (RSWT) [3] or RPA [5] as compared with 
parameters derived from mapping of microscopic models (see 
sect. 3). In all LDA(GGA)+£/ shown calculations a value U — 
J — 5 eV appropriate for edge-shared cuprates has been used. 





Ji 


•h 


a 


■h 




BSWT [3] 


-1.6 


5.59 


3.49 


0.01 


-0.4 


RSWT [3] 


-1.6 


3.56 


2.23 






RPA [5] 


-2.4 


3.4 


1.42 






present work 


-6.95 


5.2 


0.75 






3d02p optics 


-6.31 


5.05 


0.8 






GGA+U [15] 


-2.7 


15.0 


5.55 


-1.31 


0.16 


LSDA+[/ 


-8.5 


7.05 


0.82 






GGA+U 


-6.4 


5.45 


0.85 







in terms of an inproper handling of strong quantum fluc- 
tuations (SQF). Clear evidence for SQF comes from the 
small values of the ordered magnetic moment of 0.31/ie 
and the low Neel-temperature T N as 2.4 K [2, 10]. 

2. PHENOMENOLOGICAL ANALYSIS. 

The low-energy inelastic neutron scattering (INS) data [3] 
have been fitted in terms of simple (bare) 3D spin wave 
theory (BSWT), i.e. without any renormalization or ac- 
count of quantum fluctuations (see Figs. 2 (a-c,d,e)). As 
a result one arrives at jf SWT as 5.6 meV and J^swt ^ 
1.6 meV despite weak interchain couplings among them 
-0.4 meV has been claimed to be the predom- 
inant interaction responsible for the in-phase ordering of 
spirals in the magnetically ordered state below Tn (see 
Tab. 1). Then based on these and recent high-energy data 
(20 meV> u > 0.5 meV data) [5] analyzed by means of 
a random phase approximation (RPA) approach for the 
account of the coupling between the AHC, an effective ID 
model has been proposed in Ref. [5]: 



jBSWT 



ff,2 



(2/n)J. 



SWT 
2 



and Je^i as Jf 



SWT 



2J, 



SWT 



(1) 



Thereby quantum fluctuations have been taken into ac- 
count by the prefactor 2/ir as for free AHC in accord with 
the implicite a priori assumption of a weak coupling be- 
tween the AHC, in other words any renormalization re- 
lated to the coupling J e ff,i between the AHC has been 
almost ignored. Such an assumption seems to be to not 
justified from a general many-body theory point of view. 
However, following our recent work [16] we will show here 
in more detail that the in-chain exchange integrals J\ and 
J2, respectively, arc significantly different from those sug- 
gested in Re fs. [3,5,6,15]. Furthermore it has been claimed 
that the FM coupling Ji can be fixed at its bare value 

Jf SWT = jRSWT This wouW yidd jRSWT „ 3 57 mc y 

close to the result of a phemonelogical RPA-based descrip- 
tion of the problem: J^ A «3.4 meV [5]. In passing 
through we note that the value predicted by Koo et al. 



[15] exceeds that value very much by a factor larger than 
three. If the interchain coupling is of less relevance for 
the "high-energy" physics, the claimed 2J5 as —0.8 meV 
should be added to Jf PA = — 2.4 meV only for low-energy 
problems such as thermodynamics, i.e. relevant for the sat- 
uration field and the magnetization or the determination 
of the spiral's pitch angle. With such a more convinc- 
ing empirical RPA affected renormalization one would al- 
ready arrive at a = 1.063 close to the strong coupling 
boarder line. Up to now all considerations were based on 
the assumption that the FM J\ remains fixed. However, 
field-theory flow-equations based approaches [17] valid at 
a S> 1 point to strong coupling renormalizations. As a 
consequence, J\ might change considerably and a is fur- 
ther scaled down. In fact, such a tendency would be com- 
patible with our DMRG [18] results (see also below and 
Tab. 1) : J ef r,i = -6.95 meV, J cff , 2 = 5.2 cV, and a as 0.75 
[16]. If one adopts the B SWT-parameters as a reasonable 
starting point, our results should be interpreted as a strong 
upward renormalization of both | J\ | and a moderate for J2 , 
too, 

Turning first to the low-energy INS-data [3], we start 
with the two extrema of a one magnon excitation f2_ and 
f2+, i.e. the peak positions near the transferred momenta 
k = 1/4 and k = 3/4. This is the lowest two-spinon exci- 
tation (2SE) reproduced approximately by the BWST-fit 
taken from Ref. [3,5]. Its dispersion is sketched by the 
red curves in Fig. 2 (a-c). Although the maximum cor- 
responding to f2+ is broad, the asymmetry with respect 
to f2_ quantified by the dynamical asymmetry parameter 
p = f2 + /0_ > 1 is clearly visible in the INS-data in con- 
trast with the set proposed in Ref. [15] where p as 1 would 
occur. On absolute scale a discrepancy by factor exceeding 
three between the experimental and the predicted disper- 
sion is observed (see Fig 2(c) which can be traced back to 
artificially large values of J2 (see Tab. 1). 

From the experiment, Fig. 2 (d), one reads off 0_ = 
4.84 meV. Taking f2 + = 6.4 meV one estimates p as 1.32. 
p can be obtained from fitting our dynamical DMRG [19] 
results for 0.3 < a < 3 and long chains with L = 96 sites 



- + 0.0338x - 0.302^ + 0.0831x d - 0.00699a: 4 



- - 0.143a: - 0.534a^ + 0.279a; d - 0.0589a: 4 
+0.00465a; 5 , 



(2) 



where the coupling strength x = 1/a = |Ji|/</2 has been 
introduced. The relation a = f(p) provides a convenient 
highly sensitive measure of the interaction regime which is 
heavily affected by the strong quantum fluctuations. The 
function a(p) is depicted in Fig. 2 (h). One realizes ex- 
cellent agreement with a — 0.75 derived in our previous 
paper where instead fl- and the relative magnetization 
curve M(H)/M S as a function of H/H s at low tempera- 
ture have been employed [16]. Notice the large deviations 
if the BSWT or the RSWT would be applied to extract a. 
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Taking SWT w 6.7 meV and n^ SWT « 4.75 meV from 
Figs. 2 (d,e) which yields o; BSwrr « 2.3 almost consistent 
with 2.2 stated in Rcf. [3]. Using the RPA-derived values 
one arrives at a RPA about 1.42 again in formal consis- 
tency with [5] . The strong deviations of both values from 
our DMRG-based value clearly show the inapplicability 
of simple spin-wave theory based estimates. The physical 
reason is the incorrect treatment of strong quantum fluc- 
tations in the title compound which manifest themselves 
also in a small magnetic moment as mentioned above and 
in relatively large pitch angles (see below). 

Finally, considering briefly the calculated and the exper- 
imental INS intensities, at present only few comparisons 
are possible due to lacking publication of experimental 
spectra. Nevertheless, comparing e.g. the available data 
shown in Figs. 2 (d,e) one realizes that our set provides 
a better description of the intensity at large transferred 
momenta (Fig. 2(b)) as compared with that of Rcf. [5]. A 
comparison of the theoretical shapes with more INS spec- 
tra would be helpful to reduce our error bars. 

If one adopts that the experimental magnetization data 
up to the so-called field H c3 w 40.5 ± 0.2 T (H || c) where 
the peak in dM/dH occurs [12] is well-described by an ef- 
fective ID model, one arrives at the curves shown in Fig. 
3. Notice the strong deviation of the weak-coupling pro- 
posal by Koo et al. [15]. Then, a dominant FM interchain 
coupling as proposed in Rcfs. [3, 5] cannot be reconcoilcd 
with these experimental data since for such a coupling 
the 3D saturation field is smaller than its ID counter- 
part [20]. Also the spin susceptibility x(T) is within an 



RPA approach for the interchain coupling best described 
by a total AFM interchain coupling. Anyhow, a detailed 
discussion of Xs(T) including also a consideration of the 
background susceptibility \o will be given elsewhere. 

As mentioned above the presence of strong quantum 
fluctuations is evidenced by the small magnetic momen- 
tum of 0.31/i B and a low Neel-temperaturc Tn = 2.4 K. 
Both values should be compared e.g. with three to four 
times large values for the sister compound Li2CuC>2 [21,22] 
caused by a relatively strong interchain coupling [23]. In 
addition its small a « 0.32 is also helpful to suppress SQF. 
As a consequence the spiral state is significantly driven to- 
wards almost decoupled AHC the corresponding collinear 
Neel state, of each AHC i.e. the experimental pitch an- 
gle 4> = 84° analyzed within the BSWT or RSWT results 
in strongly overestimated a- values. This is illustrated in 
Figs. 4 and 5, where the maximum of the static magnetic 
structure factor is depicted as a function of a for the cases 
of a single frustrated J1-J2 chain and a coupled pair of 
them, respectively. Already the latter is expected to pro- 
vide a reasonable insight into the real quasi-lD situation. 
This point of view is supported by a detailed comparison 
with coupled cluster calculations to be reported elsewhere. 
Thus, for instance in the case of a planar arrangement of 
chains (i. e. a dominant 2D-interchain coupling as in the 
model adopted in Refs. [3,5,11,12]) the effective interchain 
interactions J| and J| correspond approximately to 

J 5 *=2J 5 ; J1 = 2J 4 . (3) 

Notice the striking failure of the classical curve especially 
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Figure 3: (Color) Magnetization vs. applied magnetic field. 
Experiment (o from Ref. [12]), theory: blue line - set of Ref. 
[5] , red line - our set [16] and black line- the (7 = 5 eV- set from 
Ref. [15] and g — 2.3 for all sets. The DMRG calculations were 
performed for L = 512 sites at T = 0. Inset: 'entire' field 
range. 



Figure 4: (Color) Pitch angle for a single chain from the maxi- 
mum of the static magnetic structure factor S(q) in comparison 
with experiment, Refs. [3,5,6], and Eq. (4) (black line). 
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Figure 5: (color) Pitch angle for two coupled chains with vari- 
ous types of interchain coupling in units of | Ji| (see Fig. 1(b)). 
The corresponding propagation vector has been estimated from 
the maximum of the calculated (DMRG) static magnetic struc- 
ture factor S(q) for L — 192 sites in total. The experimental 
value amounts about 0.467 (dashed line) [2]. 
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Figure 6: (Color) Frustration parameter a vs. effective on-site 
repulsion U for two different crystal structure refinements from 
x-ray (XRD) and neutron diffraction (ND). For each structure, 
two versions of the exchange correlation potential have been 
applied (LDA+U and GGA+C7) within the around mean field 
double counting scheme and using J-sd = 1 eV. Inset: The range 
of the obtained a values from our calculations in dependence 
of U compared with the results of Koo et al. [15]. 



84° result from the weak interchain coupling ignored in Eq. 
(4) for the sake of simplicity) and 87.42° for a = 5.55 [15]. 
What matters here is not the absolute value of </>, but the 
difference tt/2 — <j> which differs by two orders of magnitude 
between the quantum and the classic case [26]. Thus, the 
attempt to describe the spin dynamics in a quasi-classical 
way is the main reason for the inproper assignment in Ref. 
[3]. 



for large a (see Fig. 4) . Such an effect was first addressed 
in Ref. [24] in the ID-case by means of DMRG and in 
Ref. [25] for a plane of perpendicularly coupled chains 
by means of a coupled cluster approach. We stress that 
the experimental value of the pitch [2] is reproduced for 
a < 1, only, independently of the details of the weak inter- 
chain coupling. Adopting the in-chain parameters and the 
leading interchain coupling J5 w —0.4 mcV suggested in 
Ref. [3], one estimates from Fig. 5 a pitch angle of 89.58° 
(a = 2.22) and of 89.43° in the case of the RPA derived 
set ( a = 1.42) in contrast to 84° known experimentally 
[2]. Thus, the measured <j> points clearly to a strong cou- 
pling regime in contrast to opposite statements in Refs. 
[3,5,15]. Naturally, the SWT derived J's obey nearly the 
classic relation, only, (i.e. completely ignoring the strong 
quantum fluctuations) 

= cos _1 (O.25/a), (4) 

yielding 83.53° for a = 2.22 [3] (the small deviations from 



3. MICROSCOPIC ANALYSIS. Turning to a 
microscopic analysis, we compare our (DMRG) INS de- 
rived J's with those from two independent microscopi- 
cal approaches: (i) analyzing high-energy spectra from 
EELS, optical conductivity cr(uj) or RIXS data within 
strongly correlated extended multiband Hubbard mod- 
els and a subsequent mapping of their spin-states onto 
the corresponding states of a spin-Hamiltonian, i.e. the 
ID J1-J2 model under consideration. The results are 
shown in Tab. 1 and Fig. 6. (ii) extracting these ex- 
change parameters from total energy calculations of var- 
ious prepared artificially magnetically ordered states (see 
e.g. [15]). A mapping from a Cu3d 02p five-band Hub- 
bard model with usual parameters which describes the In- 
dependent dielectric response [27,28] onto a J1-J2 spin- 1/2 
model yields a sizeable J\=— 6.3 meV and J2=5.05 meV. 
We stress that in all closely related sister compounds 
[29] with a Cu-O-Cu bond angle -95° sizeable FM |Ji|- 
values 3> 1.6 meV have been found in fitting various data: 
Li 2 Cu0 2 : Ji = -1 9.6 meV (INS [23]), Ca 2 Y 2 Cu 5 Oio: 
Ji = - 14.7 mcV (INS [30]), Li 2 ZrCu0 4 : J x = -23.7 meV 
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(x(T), c p [14,31]). In particular, also for L1VC11O4 the 
T-dependent optical conductivityidata [27] obtained from 
cllipsometry measurements can be well fitted within a five- 
band Cu 3d O 2p extended Hubbard model on chain- 
clusters with up to six Cu04-plaquettes connected by 
edge-sharing. Thereby Ud — 8 cV, U p = 4.1 cV K pd = 
65 meV, A p d = 3.82 eV etc. has been used. As a result 
one arrives at in-chain J's close to the INS-derived ones: 
a = 0.8 and J2=5.1 meV (see Tab. 1). The value of J\ is 
sensitive to the magnitude of the direct FM exchange K p d 
whereas J2 is mainly sensitive to the in-chain 0-0 transfer 
integrals. Thereby |Ji| oc K p d holds approximately. No- 
tice that the contribution of K p d is much more important 
for the large negative (FM) value of 3\ than that of the 
intra-atomic FM Hund's rule coupling on O. In the past 
K p d has been used mostly as a fitting parameter rang- 
ing from 50 to 110 meV for CuGc0 3 [32,33]. A reliable 
Ji-valuc derived from an INS analysis as reported here is 
helpful to restrict its value and opens a door for systematic 
studies of this important exchange and useful comparisons 
with other sister compounds [29]. In fact, our empirical 
value corresponds to 130 meV for a er-Cu-0 bond to be 
compared with 180 meV estimated for that case and a 
cuprate plane in high-T c superconductors [34] . 

Considering the total energies of various magnetic states 
the main exchange integrals can be also extracted from 
LDA+C/ or GGA+C/ calculations. Thereby the results de- 
pend mainly on a single parameter U = U^d — J3d, where 
Jzd ~ 1 eV denotes the Hund's rule coupling that is rather 
precisely known for transition metals. For both approxi- 
mations, we calculated the exchange integrals J\ and J2 
and their ratio a for the two different crystal structures 
refined from x-ray diffraction and neutron diffraction (la- 
beled XRD and NRD, respectively, in Fig. 6 and below). 
As the key parameter, the resulting a for different values 
of U is presented in Fig. 6 and given in Table 1 (U = 5 
eV) . The graph indicates only small differences for the two 
crystal structure solutions, but essentially no difference 
for the two choices of the exchange-correlation potential 
(LDA+C/ vs. GGA+C/). For realistic parameters which 
describe successfully other edge-shared chain cuprates one 
arrives again at a < 1 in sharp contrast to Koo et al. [15] 
who obtained unusually large J 2- and a-values not com- 
patible with the observed pitch angle [2], the restricted 
two-spinon continuum, and an obviously asymmetric INS 
spectrum [5] . Presumably it is a consequence of the double 
counting procedure employed in Ref. [15] and not an arti- 
fact of the GGA as stated there because our calculations 
shown in Fig. 6 yield close values in the a-region of inter- 
est, both for the LDA and the GGA. Also the RPA-derived 
value a RPA ~ 1.4 [5] could be approached for unrealistic 
small [/-values below 3 cV adopting the XRD data, only. 

4. SUMMARY. — The main result of our revisited 
analysis of LiVCu04 is the clear evidence for strong cou- 
pling of AHC as derived from four independent experi- 
mental and theoretical studies: the INS yields a dynami- 



cal asymmetry parameter p and a pitch angle very sensi- 
tive to quantum fluctuations. Weak coupling would result 
in a nearly collinear incommensurate state and in almost 
vanishing dynamical anisotropy, i.e. p — > 1 not compati- 
ble with the diffraction and INS data. The obtained values 
for the main exchange integrals are supported by inde- 
pendent microscopic calculations based on the L(S)DA+C/ 
approach and the multiband Hubbard model. 
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